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Abstract 

We give predictions for the partial wave amplitudes of pion photoproduction 
near threshold by means of dispersion relations at fixed t. The free parameters 
of this approach are determined by a fit to experimental data in the energy 
range 160 MeV< <420 MeV. The observables near threshold are found to 
be rather sensitive to the amplitudes in the resonance region, in particular to 
the A (1232) and N* (1440). We obtain a good agreement with the existing 
threshold data for both charged and neutral pion production. Our predictions 
also agree well with the results of chiral perturbation theory, except for neutral 
pion production off the neutron. 

I. INTRODUCTION 

The photoproduction of pions near threshold has been a topic of considerable experi- 
mental and theoretical activities over the past years, ever since the Saclay |I| and Mainz 

groups showed that the data for the reaction 7p — > pvr^ were at variance with the pre- 
dictions of a low energy theorem (LET), which was derived in the early 70's by de Baenst 

01 and Vainshtein and Zakharov Being based on fundamental principles like Lorentz, 
gauge and chiral symmetry, this LET predicted the threshold amplitudes in form of a power 

1 



series in = m-^jm^^ with m^r and mj^ the masses of pion and nucleon, respectively. The 
discrepancy between the theoretical expectations and the experimental data was eventually 
explained by Bernard et al. ||^ by an explicit one-loop calculation in (relativistic) ChPT. 
It was shown that the flaw of the LET was in the assumption that the amplitudes would 
be analytical functions in /i. In fact, they contain logarithmic terms which cancel in the 
final result but whose expansion leads to previously unexpected contributions to the power 
series in /i. In the following years, calculations were considerably refined by introducing the 
heavy baryon formalism to ChPT, leading to a well-defined expansion in /i, and by including 
higher orders in the chiral expansion, up to order in the s-wave multipole £'0+ and in 
the p-wave multipoles 0. To that order there appear 3 low-energy constants, which have 
to be fitted to the data or estimated by resonance saturation. While there is some doubt 
concerning the convergence of the loop expansion for the s-waves, the expansion for the 
p-waves seems to converge rather well. Moreover, two combinations of the p-waves are free 
of low energy constants, leading to new LETs to order p'^. On the experimental side, new 
precision data have been obtained by both the TAPS collaboration at Mainz |]^ and the 
Sascatoon group. These experiments come closer to threshold than previously possible 
and clearly show that the strength of the s-wave multipole is reduced by nearly a factor 
2 as compared to the old LET. The calculations of ChPT describe the data quite nicely. 
However, there remains the question of the convergence of the loop expansion. 

Concerning the production of charged pions, the theoretical description is in much bet- 
ter shape. The dominant contribution near threshold is the Kroll-Ruderman term 0, and 
corrections up to 0(/i^) have recently been calculated in ChPT |jlO[. In that case loop cor- 



rections are fairly small, which offers the possibility to determine the pion-nucleon coupling 
constant g^^j^! by new precision experiments. 

It is the aim of this letter to investigate the threshold production of neutral and charged 
pion by means of dispersion relations at fixed t. The free parameters of this approach have 
been fitted to experiments in the energy range from 160 MeV to about 420 MeV, i.e. we 
predict the threshold behaviour on the basis of the data at the higher energies. The paper 



is organized as follows. In section ||, we summarize the kinematics and the observables 
for pion photoproduction. Some information on dispersion relations at fixed t is given in 
section |ITT|, together with a brief discussion of the assumptions and parameters involved in 
our calculations. We then discuss the predictions of dispersion theory for threshold pion 
production in section Finally, we summarize our findings and present some conclusions 
in section 0. 



II. THRESHOLD PION PHOTOPRODUCTION 



We consider the reaction 



j{q)+N{p) ^7r{k) + N'{p'). 



(1) 



The 3-momenta of photon and pion in the cm frame will be denoted by q and k, respectively. 
We will use 6cm = to describe the scattering angle, W for the cm energy of the nN system, 
and Ey for the lab energy of the incident photon. In the following we will give the multipole 
expansion for the pertinent observables taking into account only the s-wave multipoles i?o+ 
and the 3 p-waves Ei+, Mi+ and Mi_. We note, however, that our numerical results in 



section |I^ include the higher partial waves as well. For convenience we also introduce the 
combinations 



Pi = 3^1+ + Mi+ - Mi_, Pa = 3^1+ -Mi+ + Mi_, Pg = 2Mi+ + Mi_ 



(2) 



With these definitions and approximations, the differential cross section in the cm frame 



reads [|TT| 



where 



dcr I k 



■{A + P cos 6* + C cos^ 



d^l I q 



A = \Eo+ P +i I P2 P +^ I P3 |^ 
B = 2Re(Po+A*), 

c=iPip--iP2p--iP3r, 
111 2 ' ' 2 ' ' 



(3) 



(4) 



leading to the total cross section 

(Ttot = 47r||:| (I Eo+ |' + | Mi_ +6 | Ei+ \^ +2 | Mi+ j^) . (5) 
The single-polarization observables are the beam asymmetry, 

E = rsin^(|P3r-|P2|'), (6) 

the target asymmetry 

T = 2r Im {{Eo+ + cos ePi){Ps - P2)*) , (7) 
and the recoil polarization 

P = 2rim {{Eo+ + cos^Pi)(P3 + P2)*) , (8) 

where 

\k\sme fda\-' 
2\q\ [dnj ■ ^' 

The isospin decomposition of the physical amplitudes M.i{Nt:) is given by the linear 
combinations 

M^ip.-) = V2 {Mr' - M\-') = V2 (m^ - Im\'-' + Im]^' 

(10) 



III. DISPERSION RELATIONS AT FIXED t 



Assuming analyticity and unitarity, the invariant amplitudes A[s, t) for pion photopro- 
duction may be written as the sum of pole contributions and dispersion integrals. Since 
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the Mandelstam variable t is kept fixed, the integral runs over the energy variable s from 
ttN threshold to infinity. Though it is generally difficult to prove the validity of dispersion 
relations, Oehme and Taylor |12[ have given such a proof in this case, at least for sufficiently 



small values of t. As usual in this context, the pole terms are given by the Born terms 
in pseudoscalar pion-nucleon coupling. This may seem to be at variance with the well- 
established fact that pseudovector coupling should be preferred because of chiral symmetry. 
However, both approaches lead to the same contributions at the pole, and the violation of 
chiral symmetry by using the pseudoscalar coupling for the Born terms will be removed by 
appropriate contributions of the dispersion integral. 

Starting from these fixed-t dispersion relations for the invariant amplitudes of pion pho- 
toproduction, the projection of the multipole amplitudes leads to a well known system of 
integral equations |[13[| , 



1 ^oo 

ReMiiW) = M^{W) + -Y^V Ka,{W,W')lmMiiW')dW' , (11) 



(P 



where M.i stands for any of the multipoles Ei±, Mi±, and Aif for the corresponding (nucleon) 
pole term. The kernels Ku/ are known, and the real and imaginary parts of the amplitudes 
are related by unitarity. In the energy region below two-pion threshold, unitarity is expressed 



by the final state theorem of Watson |l4 



Mi{W) =1 MiiW) I e^('5'W+"'^), (12) 

where is the corresponding ttN phase shift and n an integer. 

For completeness we note that Eq. (|12D assumes that the Compton phase may be ne- 
glected, which could introduce systematical errors at the level of one percent. The further 
numerical solution of Eqs. (|1^) with the constraints imposed by unitarity follows the method 
of Schwela et al. |jl3|]. In addition we have taken account of the coupling to higher states 
neglected in that earlier reference. At the energies above two-pion threshold up to = 2 
GeV, Eq. (^) has been replaced by an ansatz based on unitarity ||T3[. Though this ansatz 



is by no means unique, it was motivated by a comparison with the existing data at the 



higher energies. Furthermore, we assume that the contribution of the dispersive integrals 



from 2 GeV to infinity is largely dominated by t-channel exchange. In earlier references |T5 
the high-energy behavior has been simulated by contributions of Regge trajectories. This 
involved a large number of free parameters which had to be fitted to the existing data with 
large error bars. Since we will restrict our calculations to energies < 450 MeV, the main 
t-channel contribution should be due to (a fraction of) p- and cj-exchange, described by four 
coupling parameters. Furthermore, we have to allow for the addition of solutions of the 
homogeneous equations to the coupled system of Eq. ([TlD . The whole procedure introduces 
10 free parameters, which have to be determined by a fit to the data. In our data base 
we have included the recent MAMI experiments for tt° and 7r+ production off the proton 
in the energy range from 160 MeV to 420 MeV [^|16|,|l7|, both older and more recent data 



from Bonn for 7r+ production off the proton |lT8|-pO|, and older Frascati and more recent 
TRIUMF data p2[ on vr" production off the neutron. 

In cases of weak coupling between the respective channels we neglected some of the 
integral kernels Ku/ in our fitting procedure. Therefore we iterated the full integral equations 
(prp to ensure consistency of the method. At low energies the only cases where we found 
a significant discrepancy between our original solution and the result after iteration, were 
the i?o+ amplitudes for neutral pion production. For this there are several reasons. First 
of all we expect strong effects from isospin breaking and unitarity for the s-wave amplitude 
near threshold. In our approach the assumption of isospin symmetry is quite essential when 



we impose the phases of the partial waves according to Eq. (|T2D. However, We introduce 
some isospin breaking by hand via the mass difference of charged and neutral pions. Second, 
the Eq+ amplitudes of neutral pion production are sums of large contributions with strong 
cancelations (see Tab. D between the large (pseudoscalar!) Born term and the dispersive 
integrals. As a consequence slight variations of these contributions can strongly affect the 
final results. 
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IV. RESULTS 



As can be seen from Tab. |, the s-waves of charged pion production at threshold are 
mainly given by the pole term contributions, while there are only small corrections from 
the dispersion integrals. In particular there is a strong cancelation between the (large) P33 
contributions to the I = \ and I = \ components entering into the linear combination Eq_^ 
(see Eqs. (Pl]|)). In this contribution had been determined to be 0.05xl0^'^/m7r, as- 
suming only the Mi+(|) amplitude to contribute. According to our result this contribution 
is 0.09x More interesting is the fact, that the contribution of the i?i+(|) mul- 

tipole which was neglected in is much larger (0.27x 10~^/m^). We also stress that 
the A contribution determined in this way differs even in sign from the corresponding 
value of -0.57xl0~'^/m7r which has been derived in heavy baryon chiral perturbation theory 
(HBChPT) [|^. However, such differences in individual terms are not unexpected, because 
there is no unique correspondence between the graphology of perturbation and dispersion 
theory. 

The first prediction for the threshold values of charged pion production was a LET 
stating that the s-wave multipoles are given by the KroU-Ruderman term 0, E^iji^n) = 
27.6 X 10^'^/m^+, i?o+ (tt^p) = —31.7 x 10^'^/m^+, where we assumed g-^jy = 13.4. Recently, 
corrections to these values have been calculated in HBChPT, to 0{^^). As can be seen 
from Tab. |T| there is good agreement for 7r+ production between our result and ChPT, and 
both agree nicely with the experiment (see also Fig. |l]). In the case of vr" production we find 
agreement of our result with the old LET value and the result of an old experimental analysis, 
while ChPT indicates a slightly stronger Eqj^ multipole (see Tab. ||). The preliminary 
analysis of the angular distribution in a recent experiment from TRIUMF (experiment E643), 
in which the inverse reaction vr^p '-/n has been studied, resulted in the much larger value 
(—34.6 ±1.0) X 10^^/m^+ |Q. At variance with that finding, the threshold extrapolation of 
the total cross sections determined from the same experiment led to (— 32.8±0.7) x 10^'^/m^+, 
reasonably close to the theoretical expectations. As has been shown in Fig. |1|, our analysis 
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describes the angular distribution at two energies quite well, although our solution for the 



Eq^ amplitude is significantly smaller than the values according to the analysis of Ref. |^ 
We note that these threshold data were not included in our fit, and conclude that there 
must be some inconsistency in the data analysis, probably due to an underestimate of the 
statistical errors in the angular distributions. There is special theoretical interest in this 
amplitude because it allows for an independent determination of the difference ai — as of the 
pion-nucleon scattering lengths via the Panofsky ratio, P = aiji^p TT^n)/cr{TT^p 7n). 



This ratio is well determined by experiment, P = 1.543 ± 0.008 [^, and related with the 
scattering lengths by time reversal, 

(ai - = 9f I El'^^in-p) \' P, (13) 







where q =\ q \ and ko =\ /c^o | are the cm momenta of photon and neutral pion at pvr" 
threshold. Using 9q/kQ = 41.3, our value of the threshold amplitude, and the measured 
Panofsky ratio, we find ai—a^ = 0.253/171^^. This has to be compared with the value 0.274 ± 
0.012/771,^ resulting from a partial wave analysis of pion-nucleon scattering (solution 
KH80). Recently, ai — 03 has also been determined by studying the level spacing of pionic 
atoms, with a preliminary value of 0.288/m7r . 



As has been previously mentioned, an exact prediction of the s-wave multipoles for 
neutral pion production is difficult. The experimental values show a pronounced cusp effect 
at the vr"*" threshold, and the expansion in fi in HBChPT seems to converge only slowly. 
In the framework of dispersion relations, we observe a delicate cancelation between the 
large negative value of the (pseudoscalar) Born term and the contributions of the dispersion 
integrals (see line 3 of Tab.l). Our final threshold value of —1.22 (here and in the following 
in units of 10^^ /mj^) is obtained after iterating the integral equations, i.e. by inserting the 
imaginary parts of our best fit into the integrals on the rhs of Eqs. ([Tl|). As has been 
mentioned, the solutions of the best fit have been obtained by neglecting some of the weak 
couplings in the system of coupled equations. In general, the iterated solutions agreed with 
the original ones at the per cent level, thus demonstrating the validity of our approximations. 
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The only exception to this is the i?o+ amphtude for neutral pion production, where we obtain 
about —0.9 at threshold without the iteration ||2^. As is shown in Fig. 0, the iterated solution 
describes the energy dependence of the data very well, in particular the steep rise of Im£^o+ 
at 7r+ threshold gives rise to the observed cusp effect. 

In the case of the reaction 7n nvr", we predict a similarly strong cancelation of Born 
terms (—5.2) and resonance contributions, leading to a final result of +1.19 (see Tab. |). 
However, our result is considerably lower than the value of 2.1 predicted by ChPT. Tab. | 
also shows the influence of the Roper multipole Mi_ on the dispersion integrals. Its size 
and even the change of sign for pn^ vs. nvr^ is of great importance for our final result. As 
the present data for pion production from the neutron suffer considerably from systematical 
and statistical errors, the predictive power of our calculation is of course much weaker in 
this case than for the proton. 

In contrast to the s-wave amplitudes, the p-waves seem to converge much faster in 
HBChPT [Q. In particular, it has been possible to derive LETs for the combination Pi 
and P2 of Eqs. (^. Our results for these multipoles are compared to the predictions of 
ChPT and experimental analyses in Tab. |ITT|. There is general agreement, typically within 
5%. The most problematic case is again the multipole Mi_. As can be seen from Eqs. (§|-§|), 
its influence is particularly strong in the case of polarization observables. We have demon- 
strated this in Figs. |] and | by arbitrarily reducing the strength of this multipole by 50%, 
which changes even the sign of the beam asymmetry S. 

V. CONCLUSION 

In conclusion we find a good overall description of the existing threshold data, showing 
the internal consistency of the data at different energies. However, there also remain some 
open questions, particularly regarding the size of the multipole Mi- connected with the 
Roper resonance, the phases of the photoproduction multipoles at energies above the second 
resonance region, and the high-energy behaviour of the amplitudes. There is also some 
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doubt concerning the accuracy of the neutron data, and it is quite clear that the present 
predictions of dispersion theory cannot be better than the data serving as input. The present 
uncertainties concerning the size of the multipole Mi_ should be removed by the analysis of 
the recently measured beam asymmetry S for neutral pion production off the proton . 
There is also hope that the proposed experiment d{e, e'd)7r^ will shed some light on the 
threshold amplitude for neutral pion photoproduction off the neutron. This reaction will be 
an extremely sensitive test of model predictions and of the size of isospin symmetry breaking 
in threshold pion production. 
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TABLES 

TABLE L Decomposition of the Eq^ amplitudes in units of 10^^/m^+. We give the contribu- 
tions of the pole terms and of the dispersion integrals (|l^) over the respective multipoles. 







pole 




^1+ 


Eo+ 




others 


sum 




_(n7r+) 


26.84 


0.13 


0.40 


-0.21 


0.91 


-0.065 


27.99 


Eq^ 


-{pTT-) 


-30.43 


-0.13 


-0.40 


-1.84 


0.86 


0.26 


-31.67 


Eq^ 


_(pvrO) 


-7.63 


4.15 


-0.41 


2.32 


0.29 


0.068 


-1.22 


Eq^ 




-5.23 


4.15 


-0.41 


3.68 


-0.93 


-0.05 


1.19 



TABLE II. The £"0+ amplitudes for charged pion photoproduction at threshold in units of 
10^'^/m^+. Our values are compared with results from chiral perturbation theory [^] and 
data analysis Q|. We also give our predictions for the reduced p-wave multipoles in units of 
I fc II ^1 10"^/m^+. 



7P vr^n 


7?! TT p 


Eo+ m.i_ ei+ mi+ 


Eo+ mi_ ei+ mi+ 


disp. 28.0 6.1 4.9 -9.6 
Ref. |10| 28.2±0.6 
Ref. |32| 28.3±0.04 


-31.7 -8.3 -4.9 11.2 
-32.7±0.6 
-31.8±0.20 
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TABLE III. s- and p~wave multipoles of neutral pion photoproduction at threshold in compar- 
ison with predictions from ChPT (6,^| and results of data analysis The unit of the s-waves 
is 10~^ /m^+, the unit of the reduced p~waves is \k\\q \ 10-^/m^+. 







disp. 


ChPT, Refs. |,31 




Exp., Ref. 


Exp., Ref. [§ 






-1.22 


-1.16 


-1.31±0.08 


-1.32±0.11 


mi_(p7r^) 




-3.91 


-3.21 




-3.38±0.26 


ei+(p7r°) 




-0.15 


-0.11 




-0.67±0.15 


mi+(p7r*') 




7.07 


7.45 




7.44±0.04 


(mi+ — mi_ 




10.99 


10.65 




10.82±0.26 


pi(p7r°) 




10.52 


10.33 


10.02±0.15 


9.2±0.3 






1.19 


2.13 






pi{mT^) 




7.77 


7.40 










-8.77 


-8.36 
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FIG. 1. Differential cross sections for charged pion production. The data for 'jn tt p are 
taken from Ref. |24], the data for jp — > vr^n from Ref. [p^ ]. 
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FIG. 2. Real part of the E'o+C^^p) amplitude before (dotted line) and after (solid line) iteration 
of Eqs. (pAj). The dashed line is the imaginary part of this amplitude. Data from Ref. Q (circles) 
and Ref. P] (diamonds). 
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FIG. 3. Differential cross section and single polarization observables for the reaction 7p 
at = 151.69 MeV (solid lines). The dotted lines show our result after reducing the Mi 
amplitude by 50%. Experimental data from Ref. M. 
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FIG. 4. Differential cross section and single polarization observables for the reaction 771 — > vr 
at E.y = 151.69 MeV (solid lines). The dotted lines show our result after reducing the Mi 
amplitude by 50%. 
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